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Abstract 

The  purpose  of  this  paper  is  to  report  on  methods  which 
have  been  developed  or  which  are  under  development 
at  Alenia  Aeronautica  for  multidisciplinary  optimum  de- 
sign, with  particulary  emphasis  on  aerodynamic  shape 
design.  Results  of  transonic  2D  and  3D  optimisation 
problems  are  presented. 

1 INTRODUCTION 

Multidisciplinary  Design  Optimisation  (MDO)  will  pro- 
vide designers  a new  array  tools  and  approaches  that 
will  take  them  closer  to  that  elusive  goal,  an  ’’optimum 
airplane”.  MDO  has  been  defined  as  ”a  formal  design 
methodology  based  on  integration  of  disciplinary  analysis 
and  sensitivity  analysis,  optimisation  and  artificial  intel- 
ligence, applicable  at  all  stages  of  the  multidisciplinary 
design  of  aerospace  systems” . What  this  really  means  is 
that  MDO  is  a way  of  getting  engineers  from  various  dis- 
ciplines, such  as  aerodynamics,  structures,  weights,  con- 
trol systems,  propulsion,  etc.,  to  work  together.  Math- 
ematical tools,  such  as  sensitivity  analysis,  modelling 
methods,  and  optimisation  solvers,  provide  a mechanism 
by  which  this  working  together  can  be  accomplished.  The 
result  is  a process  that  can  both  reduce  the  design  cost 
and  flow- time,  and  improve  the  quality  of  aerospace  sys- 
tems. 

Since  the  early  nineties,  Alenia  Aeronautica  has  im- 
proved its  optimisation  design  capability  by  participat- 
ing to  several  ECC  funded  initiatives  on  optimisation 
(EUROPT,  EC  ARP,  MDO).  The  result  is  a highly  au- 
tomatised system  that  allows  to  solve  near  real  life  design 
problems.  This  paper  describes  in  some  details  the  op- 
timisation system  and  its  potentialities.  In  particular, 
the  multidisciplinary  (aerodynamics+structure)  optimi- 
sation design  of  a wing-body  configuration  is  presented. 
Although  the  system  we  have  developed  is  general,  em- 
phasis will  be  done  in  the  following  to  wing  like  geome- 
tries which  are  the  subject  of  the  section  dedicated  to 
applications. 

2 OPTIMISATION  SYSTEM 

From  a mathematical  point  of  view,  in  order  to  solve  an 
optimisation  problem,  the  values  assigned  to  the  design 
variables  X must  be  found  so  as  to  minimize  the  objec- 
tive function  F(X)  while  maintaining  that  the  possible 
constraint  functions  Gj(X)  are  < 0.  The  designer  must 
derive  these  functions  and  choose  the  design  variables 
that  govern  the  transformation  of  the  geometry. 


The  optimisation  framework  which  was  developed  is 
formed  by  an  optimisation  module,  an  analysis  module 
and  an  interface  module  that  handles  the  geometry  mod- 
ifications (Fig.  1).  The  interface  module  is  integrated 
within  a Multi-Model  Generator  which  provides  the  CFD 
system  with  a discretisation  of  the  wing  geometry  and 
the  CSM  system  with  a finite  element  model  of  the  wing 
torsion  box  and  the  forces  acting  on  the  model  reference 
axis. 


Figure  1:  Optimisation  system. 


The  optimisation  module  is  based  on  optimisation  rou- 
tines coming  from  the  commercial  package  ADS[l],  For 
unconstrained  problems,  a quasi-Newton  method  like  the 
BGFS  variable  metric  method  is  usually  used;  while  for 
constrained  problems,  a method  based  on  Zoutendijk’s 
feasible  direction  algorithm[2]  was  selected.  The  opti- 
misation routines  are  coupled  with  analysis  codes  which 
give  the  value  of  the  derivatives  of  the  objective  function 
and  of  the  constraints.  Based  on  these  derivatives,  the 
optimisation  routine  chooses  the  most  suitable  modifica- 
tion direction.  Aerodynamic  coefficients  axe  computed 
by  using  in-house  CFD  solvers,  while  wing  mass  is  es- 
timated through  MSC  NASTRAN  structural  optimisa- 
tion module.  The  design  variables  can  be  any  of  the 
wing  planform  defining  parameters  or  those  that  govern 
the  modification  of  selected  wing  section  shapes.  The 
parametrisation  technique  is  based  on  a shape  perturba- 
tion method. 


Paper  presented  at  the  RTO  AVT  Symposium  on  “ Aerodynamic  Design  and  Optimisation  of  Flight  Vehicles  in  a 
Concurrent  Multi-Disciplinary  Environment  ”,  held  in  Ottawa,  Canada,  18-21  October  1999,  and  published  in  RTO  MP-35. 
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3 PARAMETRISATION 

The  wing  is  defined  by: 

1-  The  planform  defining  parameters  such  area(S),  as- 
pect ratio  (AR)  and  sweep  angle  (Sweep). 

2-  The  twist  (Tw)  and  thickness  (Tc)  of  key  wing  sec- 
tions. 

3-  The  shape  of  the  key  wing  sections,  which  is  gener- 
ally given  by  the  coordinates  of  a set  of  nodes. 

Each  section  is  defined  with  respect  to  a local  frame, 
centred  at  the  section  leading  edge  and  having  the  x- 
axis  coincident  with  the  section  chord.  Normalised  co- 
ordinates are  introduced  by  scaling  each  length  with  the 
chord  length,  which  means  that  x £ [0, 1], 

In  the  parametric  model,  the  shape  is  modified  using 
the  following  representation: 

x(t)  = x*(t)  = t t E [0, 1] 

n 

z(t)  = z*(t)  + ^2  Azi  BiAt) 

j- 1 

where  (x*>z*)  represents  the  position  vector  of  the  orig- 
inal shape  and  Bjin  are  general  shape  modification  func- 
tions. The  design  variables  are  the  parameters  A zj.  The 
original  shape  is  recovered  if  they  are  all  zero. 
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Figure  2:  B-spline  functions. 

The  present  method  greatly  reduces  the  number  of 
design  variables  with  respect  to  a spline  defining  curve 
approach  and  may  lead  to  a strong  reduction  on  the 
computational  cost  of  the  optimisation.  Various  shape 
functions,  including  Hick-Henne  functions[3],  Bernstein 
polynomials  and  B-spline  functions[4],  have  been  stud- 
ied. Among  those,  B-spline  approach  is  preferred  whereas 
the  functions  have  a more  compact  support  and  allow  to 
modify  only  locally  the  shape.  Fig.  2 illustrates  the 
7th-order  uniform  cubic  B-spline  functions,  which  repre- 
sents the  default  option  implemented  into  the  optimisa- 
tion system.  Only  ten  design  variables  are  required  to 
modify  both  upper  and  lower  active  wing  section  sur- 
faces. Those  sections  without  a design  variables  specified 
value  are  modified  from  the  design  variables  associated 
to  neighbouring  sections  using  a linear  relationship. 

4 ANALYSIS  TOOLS 

4.1  Aerodynamic  analysis 

The  aerodynamic  analysis  is  performed  using  simulation 
systems  based  on  in-house  developed  Euler  and  Navier- 


Stokes  multiblock  structured  or  unstructured  solvers. 
Both  take  advantage  of  explicit  time  marching  schemes 
employing  finite- volume-based  central  difference  spatial 
discretisations.  Multiblock  solvers  are  cell-centred  based, 
while  unstructured  ones  use  a node-centred/cell- vertex 
approach.  Nonlinear  second-order  and  linear  fourth- 
order  damping  terms  are  added  for  stability  and  shock- 
capturing properties.  Time  integration  is  performed  us- 
ing multistage  algorithm.  Convergence  to  steady  state 
is  accelerated  with  the  aid  of  local  time  stepping,  im- 
plicit residual  smoothing  and  enthalpy  damping  (inviscid 
computations  only).  For  turbulence  modelling,  Baldwin- 
Barth  one-equation,  together  with  A;  — a;  and  k — Rt  two- 
equation  models  are  available.  The  unstructured  solvers 
can  accepted  structured  multiblock  grids  as  input,  which 
are  converted  into  unstructured  one-block  meshes. 

The  flow  solvers  are  implemented,  together  with  a 
mesh  generator  and  a procedure  for  domain  modelling. 
Those  allow  automatic  remeshing  of  the  new  geometry 
and  of  the  associated  field  grid,  during  the  optimisation, 
as  the  design  parameters  are  modified.  As  an  alternative, 
a Laplacian  smoothing  based  node  movement  algorithm 
was  developed.  It  has  the  capability  to  deform  a grid 
while  maintaining  the  characteristics  of  the  initial  one. 
This  allows  to  introduce  an  automatic  mesh  deformation 
tool  within  the  optimisation  process,  which  we  consider 
of  paramount  importance  in  order  to  deal  with  very  com- 
plex geometries  and  with  unstructured  grids. 

4.2  Mesh  deformation  algorithm 

In  the  following,  x*  will  indicate  the  position  vector  of 
node  i.  We  will  indicate  with  Si  the  set  of  elements  be- 
longing to  the  patch  Vi  of  elements  surrounding  a given 
internal  node  i.  fci  is  instead  the  set  formed  by  the  nodes 
on  the  boundary  of  V% . 

The  method  is  based  on  an  edge  based  data  structure 
which  has  originally  been  developed  for  unstructured  fi- 
nite volume  CFD  codes.  To  move  the  mesh,  it  is  first 
assumed  that  each  node  i of  the  grid  is  connected  to  each 
adjacent  node  j by  a fictitious  spring  under  the  force  F y 
defined  by 

F ij  = Kij  (xj  — Xt)  = Kij  Axij  (l) 

where  Kij  is  the  spring  constant  which,  in  general,  will 
depend  on  some  local  grid  features  (grid  metrics). 

The  resulting  mesh  is  the  solution  of  the  equilibrium 
system  for  each  vertex  i: 

p«  = 0 - (2) 

jefCi 

which  is  also  equivalent  to  minimize  the  energy  of  the 
overall  spring  system 

2 Kij  Axy . (3) 

In  order  to  modify  the  grid  while  maintaining  it  close 
to  the  initial  mesh  characteristics,  additional  terms  are 
introduced  which  modify  the  spring  energy  relation  as 
follows 

2 y!  Kij  (A Xij  — QijAxy)2  (4) 

J'€/Ci 
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where  x°  is  the  vector  position  of  node  i in  the  initial 
grid  and  Q is  a mesh  adaptation  transformation  matrix, 
which  is  related  to  the  change  between  the  original  and 
deformed  surface  shapes.  In  the  case  of 

Q = I and  A Xij  = A x-  , 

the  spring  energy  is  zero  and  the  grid  nodes  do  not  move. 

A conjugate  gradient  method  is  used  to  minimize  the 
modified  spring  energy  equation  4.  The  technique  is  valid 
for  both  regular  and  irregular  grids. 

Fig.  3 illustrates  the  application  of  the  node  movement 
method  to  the  grid  deformation  which  results  from  a 50% 
thickness  increase  of  the  NACA0012  airfoil. 


The  last  inequality  is  introduced  to  ensure  that  the  search 
vector  remains  bounded.  Solving  this  problem  gives  a 
search  direction  which  is  tangent  to  the  critical  constraint 
boundaries,  unless  the  objective  can  be  reduced  more 
rapidly  by  moving  away  from  one  or  more  constraints. 

In  the  optimization  process,  the  gradients  are  calculated 
by  a first  order  forward  finite  difference  unless  a variable 
is  at  its  upper  bound.  In  this  case,  a first  order  backward 
finite  difference  step  is  used. 

6 APPLICATIONS 

6.1  Single-point  airfoil  optimisation 
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Figure  3:  Mesh  deformation:  50%  increase  in  thickness. 


4,3  Structural  analysis 

Starting  from  the  finite  element  and  loads  models  pro- 
vided by  the  Multi-Model  Generator,  the  wing  mass  is 
computed  using  the  structural  optimisation  module  of 
the  MSC/NASTRAN  code.  The  Sequential  Quadratic 
Programming  minimisation  method  is  usually  selected. 


5 OPTIMISATION  ALGORITHM 


At  each  iteration  q of  the  minimization  process,  the 
design  vector  X is  updated  according  to  the  formula 
X9  = X9_1-|-cjS9  , where  Sq  is  a unit  vector  representing 
a search  direction  in  a space  having  as  many  dimensions 
as  the  design  variables,  and  where  u)  defines  the  displace- 
ment in  the  direction  Sq . In  the  Method  of  Feasible  Di- 
rection, the  search  direction  is  built  in  such  a way  that  it 
will  reduce  the  objective  function  without  violating  the 
constraint  for  some  finite  move,  and  thus  can  be  defined 
as  a constrained  steepest  descent  direction.  The  search 
direction  S is  determinated  by  considering  the  following 
problem: 

Minimize: 

VF(X)  ■ S 

Subject  to: 


VGj*(X)  • S < 0 
SS  < 1 


The  TE3  test  case  [5]  adresses  the  problem  of  drag  min- 
imisation in  transonic  flows.  Far  field  conditions  are 
Moo  = 0.73  and  an  angle  of  attack  a = 2°.  The  RAE2822 
airfoil  has  been  used  as  the  starting  profile.  The  objec- 
tive was  to  search  for  a minimum  drag  based  on  a given 
lift  coefficient.  Although  this  test  case  was  defined  as  an 
inviscid  one,  optimisation  on  the  basis  of  a viscous  flow 
has  been  additionally  performed.  Two  different  test  cases 
were  run: 

1.  An  inviscid  optimisation  based  on  the  above  flow 
parameters  and  the  constraint  that  the  lift  coefficient 
for  a = 2°  ( ci  = 0.865)  had  to  be  kept  unchanged. 

2.  A viscous  optimisation  with  the  same  flow  conditions 
as  above,  a lift  coefficient  of  ci  = 0.648,  a Reynolds 
number  Re  = 6.5 -106  and  a transition  fixed  at  x/c  = 
0.03. 

Additional  constraints  have  been  applied  to  both  test 
cases.  These  constraints  read:  t/c  = (£/c)rae2822,  ol  is 
kept  fixed. 


timised  Contours. 

Inviscid  Optimisation:  The  computation  has  been 
performed  using  an  unstructured  grid  of  7809  triangles 
and  4009  nodes.  Instead  of  rebuilding  the  grid  each 
time  the  airfoil  shape  is  modified,  mesh  deformation 
was  provided  by  the  Laplacian  smoothing  based  node 
movement  algorithm.  The  ’’best  design”  was  reached 
after  69  functional  evaluations  resulting  in  a reduc- 
tion of  the  initial  drag  coefficient  which  was  driven 
down  from  its  initial  value  of  0.00746  to  0.00096.  The 
new  shape  shows  a redesign  of  the  camberline  with 
very  high  negative  value  of  the  camber  and  a thinner 
thickness  distribution  in  the  first  40  % of  the  chord, 
while  the  thickness  increases  in  the  remaining  60  %. 
The  maximum  thickness  location  moves  towards  the 
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trailing  edge  to  reach  nearly  40  % of  the  chord  (38 
% for  RAE2822  airfoil).  The  Cp  distribution  on  the 
leeward  side  of  the  optimised  airfoil  is  characterised  by 
an  increase  of  the  suction  peack  and  a reduction  of  the 
shock  jump;  the  winward  side  remains  almost  unchanged. 


Figure  5:  Problem  TE3  - Transonic  Inviscid  Drag  Min- 
imisation - Comparison  of  Pressure  Distributions. 

Viscous  Optimisation:  The  computation  has  been 

performed  on  an  almost  structured  grid  of  8200  elements 
and  8360  nodes.  Here  again,  mesh  deformations  were  ob- 
tained using  the  node  movement  algorithm.  The  ’’best 
design”  was  reached  after  68  functional  evaluations  re- 
sulting in  a reduction  of  7 drag  counts  ( Cd  = 0.0081  in- 
stead of  0.0088).  The  shape  is  closer  to  the  RAE2822  air- 
foil than  those  coming  from  inviscid  optimisation.  Here 
again,  reduction  of  drag  is  mainly  obtained  by  the  re- 
design of  the  camberline  but  in  a lower  measure.  The 
thickness  distribution  remains  almost  unchanged.  The 
Cp  distribution  is  characterised  by  a slightly  higher  suc- 
tion peack  and  is  shock  free. 


Figure  6:  Problem  TE3  - Transonic  Viscous  Drag  Min- 
imisation - Comparison  of  Pressure  Distributions. 

6.2  Multi-point  airfoil  optimisation 

The  problems  consist  in  the  improvement  of  an  already 
existing  airfoil  in  order  to  increase  its  maximum  speed 
without  deteriorating  its  characteristics  at  lower  speeds. 
More  precisely,  the  design  problem  was  defined  as  to  min- 
imize drag  at  Moo  = 0.8  and  Ci  = Cf , maintaining  the 
Ci  and  Cd  characteristics  on  the  buffet  onset  curve  of  the 
original  airfoil  for  Moo  = 0.55,0.62,0.71,0.754,  which  re- 
sults in  a five-point  optimisation  problem.  The  value  of 


the  Mach  number  is  24  • 106  and  the  flow  is  assumed  to 
be  fully  turbulent. 

The  new  shape  is  obtained  by  a redesign  of  the  camber- 
line with  higher  values  of  the  camber  on  both  leading  and 
rear  parts  of  the  airfoil,  while  the  thickness  is  increased 
beyond  40%  of  the  chord.  The  maximum  thickness  posi- 
tion moves  from  35%  to  42%  of  the  airfoil  chord. 
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Figure  7:  Multi-Point  Optimisation  - Comparison  of  Ini- 
tial and  Optimised  Contours. 
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Polar  at  M =0.55. 


Polar  at  M=0.62. 

Figure  8:  Multi-Point  Optimisation  - Comparison  of 

Aerodynamic  Characteristics. 

Figures  8 and  9 illustrate  the  aerodynamic  characteris- 
tics ( Ci  — Cd  polars)  of  the  two  airfoils.  At  the  main  de- 
sign point  (M  = 0.80),  the  drag  is  reduced  along  the  over- 
all lift  range.  For  the  secondary  points,  a small  penalty 
on  the  drag  (a  few  counts)  is  observed  for  low  and  in- 
termediate values  of  the  lift,  while  the  drag  is  reduced 
at  higer  values.  In  general,  an  increase  of  the  maximum 
lift  of  the  airfoil  is  reached.  This  is  particularily  true  at 
M = 0.71  and  M = 0.754. 
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Polar  at  M=0.80. 


Figure  9:  Multi-Point  Optimisation  - Comparison  of 

Aerodynamic  Characteristics. 

6.3  Transonic  wing  optimisation 

The  reference  aircraft  geometry  selected  within  the  MDO 
project  [6]  is  representative  of  a wide-body  civil  transport 
aircraft  (650  passengers),  with  a wing  span  of  nearly  80 
meters  and  a maximum  take-off  weight  of  about  550  tons. 
The  wing,  whose  exposed  right  component  planform  is 
illustrated  in  Fig.  10,  has  an  area  of  725  m2,  an  aspect- 
ratio  of  8.2,  a span  of  77.1  m and  a sweep  angle  of  33°. 


The  optimisation  problem,  which  has  been  solved,  is 
summarized  as  follows: 


Design  Operation 

Economic  Cruise(Moo  = 0.85, 
Cfot  = 0.458) 

Objective  Function 

ADOC  = AD  + AW 

Constraints 

Cl  = cf 
Cm=Ctm 

Design  Space 

Crank  and  Tip  section  shapes; 
Crank  thickness 

where  AD  is  the  change  of  the  wing-body  drag  expressed 


in  counts  and  AW  is  the  change  in  mass  of  one  torsion 
box  expressed  in  tons. 

Starting  from  a reference  trimmed  configuration,  the 
problem  of  reducing  the  DOC  by  only  operating  on  the 
wing  geometry  without  modifying  other  parts  of  the  air- 
plane is  thus  addressed.  The  trim  condition  is  satisfied 
by  maintaining  the  same  angle  of  attack,  the  same  lift 
and  the  same  pitch  moment  than  those  of  the  original 
configuration.  The  root  section  is  kept  fixed.  Crank  and 
tip  wing  sections  just  as  the  crank  thickness  are  allowed 
to  change,  which  corresponds  to  21  design  variables.  In 
addition,  it  is  assumed  that  crank  thickness  modification 
induces  the  same  percentage  change  in  tip  thickness. 

The  structural  wing  model  consists  in  a centre  wing 
box  which  is  integrally  built  into  the  centre  fuselage  to 
which  are  attached  the  port  and  starboard  lateral  wing 
sections.  The  centre  wing  box  is  a conventional  metallic 
structure  with  three  spars.  The  structural  optimisation 
model  of  the  wing  uses  156  design  variables  involving 
size  parameters.  The  finite  element  model  is  formed  by 
around  12000  elements.  A fine  grid  of  54  blocks  and 
about  204000  cells  has  been  built  around  the  wing  body 
configuration.  It  is  used  to  compute  the  trim  condition 
of  the  reference  geometry  and  to  verify  the  aerodynamic 
characteristics  of  the  optimised  one.  A coarse  grid,  which 
is  the  eighth  of  the  fine,  has  been  implemented  into  the 
optimisation  system  in  order  to  maintain  the  overall  opti- 
misation cost  at  a reasonable  level.  For  the  same  reason, 
the  optimised  wing  mass  was  approximated  by  a curve 
fit  as  function  of  the  crank  thickness. 


Table  1:  Aerodynamic  performance  of  the  optimised  con- 
figuration^  


Case 

Config. 

Cl 

Cm 

Cd 

Ref 

w-b  coarse  grid 

0.4946 

-0.1963 

0.02120 

w-b  fine  grid 

0.5131 

-0.2047 

0.01700 

exposed  wing 

0.4550 

-0.1835 

0.01390 

Opt 

w-b  coarse  grid 

0.4928 

-0.1957 

0.01790 

w-b  fine  grid 

0.5154 

-0.2049 

0.01370 

exposed  wing 

0.4579 

-0.1848 

0.01060 

The  reduction  of  drag  only  results  from  a reduction  of 
the  wing  wave  drag,  as  it  was  expected  from  the  problem 
definition.  Aerodynamic  performances  of  fuselage  and 
tail  are  maintained  unchanged.  The  use  of  a coarse  grid 
seems  to  be  adequate  since  the  drag  variation  between 
the  reference  geometry  and  the  optimised  one  is  nearly 
the  same  for  the  two  grids. 

The  improved  configuration  is  obtained  after  a few 
global  iterations  of  the  optimisation  process,  which  never- 
theless corresponds  to  around  150  functional  evaluations. 
The  DOC  is  reduced  while  meeting  the  constraints  at 
each  iterations. 
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Table  2:  Wing  section  thickness  (%)  and  wing  mass  (Kg) 
of  the  optimised  configurations. 


Case 

RootTc 

CrankTc 

TipTc 

Wing  Mass 

Ref 

14.00 

10.00 

8.800 

28764 

Opt 

14.00 

8.785 

7.730 

30562 
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Figure  11:  MDO  wing  - Comparison  of  Initial  and  Opti- 
mised Contours. 


For  the  optimised  geometry,  the  upper  surface  of  the 
crank  section  is  flattened.  This  effect  has  been  mainly 
obtained  by  the  redesign  of  the  camber  line.  Except  the 
decrease  in  the  maximum  thickness,  the  thickness  dis- 
tribution is  only  slightly  altered  in  the  rear  part  of  the 
airfoil.  In  the  case  of  the  tip  section,  the  major  changes 
again  apply  on  the  camberline,  with  a tiny  reduction  of 
the  leading  edge  curvature  in  the  thickness  distribution. 
The  sectional  pressure  plots  indicate  a movement  towards 
a two-shock  structure  on  the  upper  surface. 


ETA 


ETA 

Figure  12:  MDO  wing  - Comparison  of  Span  Loads. 
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Figure  13:  MDO  wing  - Comparison  of  Pressure  Distri- 
butions. 


crank  thickness  modification,  the  reduction  of  crank  sec- 
tion thickness  has  the  effect  of  further  reducing  drag  not 
only  in  the  central  part  of  the  wing,  but  also  inboard.  It 
leads  nevertheless  to  an  increase  of  the  wing  mass. 

7 CONCLUSION 

Mono-  and  Multi-disciplinary  capability  of  the  optimi- 
sation system  has  successfully  been  demonstrated  for 
2D  and  3D  applications.  The  capability  of  the  design 
tool  to  be  able  to  account  for  explicit  constraints  and 
multi-objectives  within  the  optimisation  process  is  of 
paramount  importance  for  industrial  applications.  Nev- 
ertheless, the  cost  of  sensitivities  is  a key  factor  when 
advanced  modelling  methods  are  used.  Cheaper  sensitiv- 
ities and  parallel  implementations  are  requested.  Future 
works  will  involve 

1.  more  accurate  gradients  computation  (adjoint  equa- 
tions, use  of  automatic  differentiation), 

2.  introduction  of  other  disciplines  within  the  design 
loop  (aeroelasticity,  electromagnetism,  . . . ). 

8 ACKNOWLEDGEMENTS 


Spanwise  loads  distributions  display  a decrease  of  the 
drag  in  the  central  part  of  the  wing,  where  lift  and  pitch 
moment  are  approximatively  maintained  at  the  same 
level  than  those  of  the  reference  wing.  The  loss  of  lift 
in  the  inboard  part  of  the  wing  is  recovered  by  the  out- 
board sections.  The  opposite  behaviour  is  found  on  the 
pitch  moment.  By  comparison  with  the  problem  without 


Part  of  this  work  was  performed  within  MDO  project: 
a collaboration  between  British  Aerospace,  Aerospa- 
tiale, DASA,  Dassault,  SAAB,  CASA,  Alenia,  Aerma- 
cchi,  HAI,  NLR,  DERA,  ONERA  and  the  Universities 
of  Delft  and  Cranfield.  The  project  was  funded  by  the 
CEC  under  the  BRITE-EURAM  initiative  (Project  Ref: 
BE95-2056). 


20-7 


References 

[1]  Vanderplaats,  G.N.,  ADS  - A Fortran  Program 
for  Automated  Design  Synthesis , User  Manual, 
1987. 

[2]  Vanderplaats,  G.N.,  Numerical  Optimization 
Techniques  for  Engineering  Design , McGraw- 
Hill  Series  in  Mechanical  Engineering. 

[3]  Hicks,  R.  and  Henne,  P.A. , Wing  Design  by  Nu- 
merical Optimization , Journal  of  Aircraft,  vol 
35,  n.  7,  1978. 

[4]  Far  in,  G.,  Curves  and  Surfaces  for  Computer 
Aided  Geometric  Design , Academic  Press,  1989. 

[5]  Periaux  J.  et  al  eds.,  Optimum  Aerodynamic  De- 
sign and  Parallel  Navier-Stokes  Computations , 
Vieweg,  1998. 

[6]  Allwright  S.,  MDO  Process  and  Specification  for 
the  Primary  Sensitivity  Study , MDO  Technical 
Report  MDO/SPEC/BAE/SA960430,  1996. 


DISCUSSION 


Session  III,  Paper  #20 


Prof  Slooff  (NLR,  Netherlands)  noted  that  in  their  work,  optimising  on  the  basis  of  Euler 
equations  often  led  to  non-unique  solutions,  that  is  several  minima  of  about  the  same  level.  In 
such  cases  the  additional  constraints  helped  to  reduce  or  eliminate  this  “non-uniqueness”.  Had 
the  authors  encountered  this  problem? 

Dr  Selmin  indicated  no  specific  experience  of  this  issue,  but  agreed  that  adding 
constraints  should  help.  He  saw  this  as  a requirement  for  a more  precise  definition  of  the 
original  design  problem. 

Mr  Perrier  (Dassault,  France)  asked  how  better  predictions  of  sensitivity  could  be  ensured  as 
he  believed  that  decisions  on  constraints  should  come  from  good  knowledge  of  constraint 
sensitivity. 

Dr  Selmin  noted  that  his  team  intended  to  improve  prediction  of  sensitivity  by  moving 
from  the  use  of  finite  differences  to  automatic  differentiation  to  evaluate  the  gradients. 


